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Abstract 



The convolution formula is derived within the framework of the decay-chain method 

m 

' for decay channels with three and four particles in a final state. To get this formula 

^3 . exactly for unstable particles of any type one must modify the propagators of vector 

00 ' and spinor fields. In this work we suggest proper modifications and get the convolution 

O '. 

, formula by direct calculations. It was noted that this approach naturally arises in the 
O 1 

-■»«. ■ model of unstable particles with random mass. 
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1. Introduction 

The unstable particles (UP) manifest themselves as resonances (intermediate states) in the 
scattering experiments and as decaying (evolving) states in the decay or oscillation experi- 
ments. The former case is treated with S - matrix or with renormalized propagator approach. 
In the second case we must take into account the instability or finite-width effects (FWE) 
in a special way. There are some different methods proposed in literature to describe the 
UP as quasi-stable state. The first method was formulated by Matthews and Salam in Ref. 
[ij, where an uncertainty principle was taken into consideration. Bohm et al. [2] suggest the 
time asymmetric quantum theory of UP with relativistic Gamov vector, which describes UP 
in initial or final states. The model of UP with random (smeared, fuzzy) mass was suggested 
in Ref. [3j. This model is based on the uncertainty relation for the mass of UP and has a 
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close analogy to the elaboration pQ. An effective theory of UP is discussed in [3] where the 
authors have applied the method of correction factorization. 

At last decade the so called convolution method (CM) has been used to evaluate FWE 
[HIE]. This method consists in a semiphenomenological description of invariant mass distri- 
bution by Breit-Wigner-like density function. The convolution formula was derived within 
the framework of the model [3] as a direct consequence of the model approach. This formula 
was applied in a phenomenological way to calculate the hadron decay rates [7] . The FWE (or 
"mass smearing" effects) in hadron decays are very significant due to large hadron widths. 
The FWE in the near-threshold decays t -> WZb, t -> cWW, cZZ and A°(h°) -> tbW were 
calculated in Refs. [3, E] where some analysis of CM applicability is fulfilled. 

In this paper we systematically analyze the processes of type $ — ► 0i0 — > 010203, 
where is an UP with a large width and $,0i,02 are the stable or long-lived particles of 
any kind. The prescriptions are suggested for propagators and for polarization matrices 
of unstable vector and spinor fields, which exactly lead to the convolution formula for a 
decay rate. It was noted that these prescriptions naturally follow from the model [3]. The 
convolution formula for above mentioned decays was got by direct calculations without any 
approximations for all type of UP and for all possible type of $ and <p a . 



2. Derivation of the convolution formula for a decay 



The convolution formula (CF) can be obtained on the example of processes $ — > 0i0 
12i k 0i0i0fc ; where is an UP with a large width. In the general case: 



where the sum runs over all decay channels of the unstable particle 0. Some particles of the 
X n = (xi,X2, ■■■,x n ) can be unstable too. By means of Eq. (DQ) we represent the connection 
between the CM and decay-chain method (DCM). When the sum runs over (0^,0^), which 
are the stable or long-lived particles, we have: 



rate 




(i) 



n 




(2) 



f k 



The CM gives the expression [5]: 




(3) 
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In Eq. (J3J) the value T($ — > 010(g)) is the width, calculated in a stable particle approximation 
when mi = q 2 and p^(q) is some invariant mass distribution function. From Eq. (PQ) it is 
clear that p^{q) depends on all decay-chain of 0, that is, CF can be derived in principle from 
the DCM. 

When the intermediate state is scalar UP then for any type of initial $ and final 
states a decay width can be represented in the factored form (see Appendix): 



hi 



m-n y( , 2 n :i )= I r(a->^( g )) 2 V|p 0(g) j 2 rf g • ( 4 ) 

In Eq. (JIJ) 9i = m 2 + m 3 , g 2 = m$ — mi and P^,(q) is the propagator's denominator of the 
scalar field 0. It is noteworthy that the derivation of the convolution formulae ([3"]) or (jlj) 
don't depends on the choice of P^>{q). This fact is important for a model definition of p(q). 

For convenience and completeness we give in Appendix a full list of the expressions 
T($ — > 0i0) and T($ — > 010203) for all possible types of $ and 0&. From Eqs. (j2J) and (jlj) 
it follows that /fy(?) in Eq. (j3J) can be expressed in the form: 

A generalization of Eq. (J5J) for the cases when three or more particles are in a final state 
is straightforward (see Appendix). Thus, we can always represent r($ — > 0i0), where is 
scalar UP with a large width, by the convolution formula ([3]). If one uses the parametrization 
qT(q) = IrriE(q) and Dyson-resummed propagator 

P M) = <? 2 - rnl(q) - UwE^q), mj(g) = + i?eS (g), (6) 

then the p<^(?) can be written in the Lorentzian (Breit-Wigner type) form: 

/ % 1 ImE^jq) , v 

^ ^ [ g 2_ m 2 ] 2 + [JmE ^ (g)] 2- U 

The expression (J7]) have been used in the previous papers El [7] . 

The situation drastically changes when UP is vector or spinor field. Traditional propaga- 
tors Pfjiv^q) = —iig^u — q^qv l m v) I Pv(q) for the vector field V and P(q) = i(q + rn^)/P^(q) 
for the spinor field \P do not make it possible to represent T($ — > 010203) in the fac- 
tored form (jlj). As it was noted in Ref. [6] for the vector field case the numerator r]^ = 
~ Q^Qf/^^v) destroys the factorization and we have an approximate CF. The same 
disfactorization takes place for the fermion propagator i(q + m^) too. In this situation 
one should analyze the modification of propagators for vector and spinor field after Dyson 
summation. Such an analysis is complicated due to tensor and operator structure of r]^ v and 



q + m. There are no unit and strict definitions of these structures in literature. For example, 
the rj^ has the structure g^ - q^q v /(w^ - im v T° v ) in Ref. [6], g^ - q^q v /(m v - TV/2) 2 
in Ref. [5], - q^q v /q 2 in Ref. [8], and (g^ - q^q u/ 'q 2 )Pr(q 2 ) + {q^u/ q 2 )PL{q 2 ) in Ref. 
[9]. By direct calculation one can check that T]^ = —i^g^u — q^qu/q 2 ) leads to the factored 
expression (JIJ), that is to CF. 

The model of UP with random mass [3] contains the general designation how to modify 
propagators P^uiq) and P(q). The "smearing" of UP's mass in according to uncertainty 
principle leads to the modified dispersion relation q 2 = fi 2 , where \i is nonfixed random mass 
of UP with some distribution function p(/j). In another way we have the "smeared" or 
"fuzzy" mass-shell and the general prescription m 2 — > fi 2 = q 2 , where q is timelike arbitrary 
momentum. Then the modified propagators for vector and spinor fields are: 



In Eqs. (JHJ) and ([9]) the Pv(q) an d P^(q) are some functions, which don't influence on 
convolution structure of the formulae ([3]) and (HI). The prescription (jHJ) coincides with the 
definition in Ref. [8]. It should be noted that the expressions (JHJ) and ([9]) approximately 
coincide with standard ones at peak vicinity q 2 ~ m 2 (m). 

Direct calculations with help of Eqs. (JHJ) and (JHJ) (see Appendix) lead to the factored form 
(JH) for r($ — > 0102^3) when the UP in an intermediate state is vector or spinor field. The 
calculations were fulfilled for all possible type of particles $ and 0&. For completeness we 
give a full list of the expressions for r($ — > 010203) in Appendix. Thus, the prescriptions (JHJ) 
and (JHJ) always lead to CF (J1J) or (J3j) and (JHJ). The expression (J7J) runs out from Eq. ([5]) and 
parametrization qT(q) = imX(g), when we use the expression ([6]). It should be noted that 
applicability of this approach in the context of gauge theories is limited to low orders because 
of the appearance of gauge dependence JHJ JHJ JTOJ JHJ EG] . Moreover, we should redefine Py(q) 
for spinor UP with account of m^(q) = m ^ + i?e£$(g) and T^(q) = JmS^(g). However, as 
was noted early, the convolution structure is not subject to definition of P(q). 

To get the expression for T($ — > 010(g)), where 0(g) is V or one need the definitions of 
polarization matrixes for vector and spinor fields. With prescription m — > q these definitions 
are: 



Piwio) = -*(<//» - q»qv/q 2 )/Pv{q) 
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and 



P(g) = j($ + q)/P iS! (q), q = y/{gqj . 



(9) 




(10) 



e 
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and 

The expressions (ISjl - (fTT]) are the continuation of standard ones to the "fuzzy" (smeared) 
mass-shell. The factorization of T($ — > 010203) on "fuzzy" mass-shell results due to the 
right sides of Eqs. (TTU]) and (TTTT) equal to the numerators of Eqs. (jSD and Such a 
factorization takes place on the usual mass-shell. 

Direct calculations do not lead to the CF in the cases when there are two or more UP in 
the intermediate or final states. In the framework of the model [3J for $ — > 0102, where 0i 
and 02 are the UP with large width, we get the expression for the width, which is usually 
applied in a phenomenological way: 

r($ - 0!0 2 ) = J dqj Pl ( qi ) J dq 2 2 p 2 (q 2 )T^ - MlMte))- (12) 

Representation of the width T($ — > 010203) in the form is actually the transition from 
DCM to CM. Such transitions are very complicated in the cases when there are many unstable 
and stable particles and there are many sections in the decay-chain. Instead of transition 
from DCM to CM we can derive the convolution formula in the framework of the model [3] 
for general case. In Ref. [3] such a derivation was done in the case of scalar UP, therefore 
the discussion of this problem for the vector and spinor UP is actual. 



3. Summary 

In this paper we have demonstrated in detail the connection between DCM and CM for the 
decay channels with three-particle final states. The convolution formula was derived with 
help of the prescription m — > q by direct calculations for all types of UP. It was noted that 
this prescription naturally arises in the model of UP with random mass [3J. Transition from 
DCM to CM is very complicated in the case of many-particle decay-chain. Therefore an 
alternative proof of CF in the general case is actual. 



4. Appendix 

In this section we use the A-function, which describes the kinematics of the process 0(g) — > 
01 (A;i)02(A;2)- In the q = frame of reference: 

ka = ^qMh, h; q),k a = \k a \, a = 1, 2 , (13) 
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where the X(k±, k 2 ; q) is defined in analogy with the known Kallen function X(k\, k 2 ; q) 
g 2 A 2 (fci, fc 2 ;g) 

1.2 1.2 l,:2 7.2x2 

g z q" 

The expressions 1^(0 — > 0i0 2 ) for all types of particles can be represented in the form: 



pip /p _ p\2 

A(*i, fc 2 ; g) = [1 - 2 + ^ , 2) ] 1/2 • (14) 



3 2 ~ 

r,(0 -> 0i0 2 ) = — A(mi, m 2 ; m)f i (m 1 ,m 2 ; m), (15) 

<37T 

where /j(rai, m 2 ; m) depends on the interaction Lagrangian. To illustrate CM we use the 
simplest Lagrangians: 

£ fc = <700i0 2 ; <?0V^ ; 0#iV> 2 ; <W M & - 0f 0); <?V^ l7 M (c y + c A76 )^; ffW^V (16) 

In Eq. ( |T6|) 0, and are the scalar, vector and spinor fields. Then the /j(m 1 ,m 2 ;m) in 
Eq. (fl5|) are defined by the following expressions: 
1 

2m 

^ VlV2i /2 = l [1+ (™ 2 -"f--"i> 2 ]. (18) 

0-« 2 , /3 = m[l- (mi+ 2 m2)2 ]. (19) 

j 777,^ ~ 

— > 0i V, U = -m — K\ 2 (m 1 ,m 2 ;m), m 2 = m v . (20) 
2 wij 

V -> 0i02, /5 = ^mA 2 (mi,m 2 ;m), m = my. (21) 

v - /. = ?m{(4 + 4)H - 4±^l _ K^ll!, + 3(4 _ (22) 



i*2, A = s-- (17) 



1 (m 2 + m? — m 2 ,) 2 
— 1 + " 

3m 8m 2 m 2 

1 m 5 



^ - ^10, A = ^[1 + 1 V 7 2 2J ], m 2 = m,. (23) 



V -> / 8 = ^-^{l + 8(/ii + /i 2 ) - 2(9/i 2 + 16/ii^ + 9/j 2 ) + 8(j4 - 4/i 2 /i 2 (24) 

24 mfm 2 
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- 4/Ji/i 2 + /i 2 ) + /4 + 8/j 2 /j 2 + 8/ii/i 2 + /i 2 }, /_f Q = — f . 



1 mi m? — m 2 
ip —> (pipt, 79 = 7:^(1 + 2 1 ^ — ), m = m^,m 2 = m^. (25) 

2 mm 1 



2 2 ( m2 — m i) 2 + 'mKrn 2 + m 2 — 2m 2 ,) 



^ -> ^iv, / 10 = (4 + 4)^ — 1 ' — ~"" 27 - (4 - 4), (26) 

m = 77ty, r«2 = my. 



Using the expressions (17) -(26) we can represent T($ — > 010203) in a compact and 
universal form for all types of decay channels. Here we shortly describe the method of 
r($ — > 010203) calculation. This value can be always written as: 

r = ^/j(|M(fc i ,m 4 )| 2 )^ I (27) 
6 



where M(ki,rrii) is an amplitude, p and ki are momentum of $ and fa, k is some numerical 
factor, and 

J(\M\ 2 ) = / \M\ 2 5(p -k x -k 2 - h)^. (28) 

The integral J(\M\ 2 ) is calculated in q = frame of reference and as a result we have the 
noncovariant expression 

J{\M\ 2 ) — f(q°,qy,q°klp 2 ,...). (29) 
This expression can be always reconstructed to covariant form using q = 0: 

q°^q= Vim), q°P° - (qp), A? - f = {p°? ~rn 2 ^ {pqf/q 2 - m 2 , ... (30) 



Then we pass to the p = frame of reference and change the variable in Eq. (1271) according 
to 

—7^- = k 1 dq 2 dfl = — \(q,mi] m)dq 2 dfl. (31) 

ki 2m 4 

Using this simple method and prescriptions ([8]), (Q we have got by tedious but straight- 
forward calculations the general expression for T($ — > 010203), when $, and 0& are all of 
possible type particles: 

T af3 ($ -> 0i0i0 fc ) = f^f / A(g, mi; m)f a (q, m x ; m)\(m h m k ; q)fp{m i , m k ; q) ^ , 

(32) 

where gi = m; + and q 2 = m — m\. From Eqs. ( 1321) and ( TT5l) it follows: 



_• 



qTp{(j>{q) fafa) 



r Q/3 ($ - 0i0 4 fc ) = / dg 2 r a ($ - 010(g)) _pr^ . (33) 

In the approximation 

r($ - 0,0) = ^ r($ -> 0i0,0 fc ) (34) 



i.k 

we get the known convolution formula 



IVI> - c _) 1. ' ; ) = r($ - 0i0(g))p^(g)rfg 2 , (35) 
9? 



where 

The same result can be received for many-particle decay channels of UP — > 010203- • • F° r 
example, let us consider the decay chain $ — > 0i0 — > 01020304, where 0& are the scalar 
fields. Then for the simplest contact interaction we have: 



9i92 f dki dk 2 dk 3 dk 



r * = 25%vi mm? J J /^ fc2 ^ 3 - M ^ppr' (37) 



where q = p — k\ and 

r*(g) = r(0( g ) -> </>i0 2 3 ) = / / Js(q-k 2 -k 3 - • (38) 

From Eqs. (S3), (EE} and (??) it follows: 

r* = [ q2 d q 2 rM f p j q ] l2 , (39) 

where T$(g) = T($ — > (pi(f)(q)). Using the factorizable |M| 2 we can get the result fl39l) by 
direct calculations for others types of particles It should be noted that the factored (1331) 
and convolution (|35j) structures take place for any choice of P,p(q)- 
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